Numerical Study of Spin Hall Transport in a Two Dimensional Hole Gas System 
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We present a numerical study of the spin Hall effect in a two-dimensional hole gas (2DHG) 
system in the presence of disorder. We find that the spin Hall conductance (SHC), extrapolated to 
the thermodynamic limit, remains finite in a wide range of disorder strengths for a closed system 
on torus. But there is no intrinsic spin Hall accumulation as induced by an external electric field 
once the disorder is turned on. The latter is examined by performing a Laughlin's Gedanken gauge 
experiment numerically with the adiabatic insertion of a flux quantum in a belt-shaped sample, in 
which the absence of level crossing is found under the disorder effect. Without disorder, on the other 
hand, energy levels do cross each other, which results in an oscillating spin-density-modulation at 
the sample boundary after the insertion of one flux quantum in the belt-shaped system. But the 
corresponding net spin transfer is only about one order of magnitude smaller than what is expected 
from the bulk SHC. These apparently contradictory results can be attributed to the violation of 
the spin conservation law in such a system. We also briefly address the dissipative Fermi surface 
contribution to spin polarization, which may be relevant to experimental measurements. 



I. INTRODUCTION 

Recently it has been proposed 1 ' 2 that in spin-orbit cou- 
pling (SOC) systems one may use an electric field to gen- 
erate transverse spin currents in the absence of external 
magnetic fields. It has been argued 1,2 that such a spin 
Hall effect (SHE) is intrinsic, contributed by all the elec- 
trons below the Fermi energy, and the corresponding spin 
currents are "dissipationless" as in contrast to the dissi- 
pative longitudinal charge currents which are only con- 
tributed by the electrons close to the Fermi energy and 
are strongly subjected to scattering effects. 

The original proposals for the SHE are for the disorder- 
free cases. 1 ' 2 In the two-dimensional (2D) electron gas 
described by the Rashba model, 2 it was shown based on 
a perturbative approach 3 ' 4 that the SHC is precisely can- 
celled by the vertex correction once the disorder is turned 
on. On the other hand, the vertex correction is found to 
vanish 5 for the three-dimensional (3D) p-doped semicon- 
ductors described by the Luttinger model 1 so that the 
SHC is still finite in the presence of weak disorder. Nu- 
merical calculations 6-8 of the SHC seem to support such 
perturbative results of the distinct behavior for two mod- 
els in the thermodynamic limit. 

Experimentally the signatures of spin polarization have 
been observed recently in 2DHG system 9 and 3D n- 
doped semiconductors, 10 which have generated a lot of 
excitement concerning whether they are due to the in- 
trinsic SHE mentioned above or some extrinsic effect. 11 
Bernevig and Zhang have shown that the vortex cor- 
rection does vanish in the 2DHG 12 and 3D n-doped 
semiconductor 13 systems. Furthermore, the mesoscopic 
SHE in the 2DHG is also found to be present based on the 
nonequilibrium Green function method 1 and Landauer- 
Biittiker formula with attached leads, 15 similar to (with 
larger magnitude than) the mesoscopic SHE found in the 
2D Rashba model. 16,17 But it is still unclear whether the 
bulk SHE in the 2DHG can survive in the thermodynamic 
limit beyond the perturbative approach. 



In this paper, we perform numerical calculations for 
the 2DHG in the presence of disorder. We first show 
that the bulk SHC calculated from the Kubo formula is 
indeed robust against the disorder in extrapolation to the 
thermodynamic limit, which is consistent with the vertex 
correction calculation. 12 It is also similar to the behavior 
for the Luttinger model, 7 but is in opposite to that of the 
2D Rashba model. 6 But when we perform a Laughlin's 
Gedanken "gauge experiment" on a belt-shaped sample 
to probe the spin transfer/accumulation due to the SHE, 
we obtain a null result due to the anticrossing between 
energy levels, which is quite similar to what has been pre- 
viously seen in the 2D Rashba model. 6 We point out that 
the absence of edge states in the 2DHG system causes the 
general level repulsion with the turn on of disorder, which 
leads to the disappearance of the net spin Hall accumu- 
lation in an open system. Furthermore, at zero disorder 
case, energy levels do cross each other and we find an 
oscillating spin-density-modulation at the sample bound- 
ary after an adiabatic insertion of a flux quantum in the 
above-mentioned Laughlin's gauge experiment. However, 
the corresponding net spin transfer is only about one or- 
der of magnitude smaller than what is expected from the 
calculated bulk SHC. We discuss the "conflicting" results 
of a finite SHC but without an intrinsic spin accumula- 
tion in the 2DHG, and point out that the underlying 
reason can be attributed to the violation of the spin con- 
servation law in such a SOC system, where the SHC is no 
longer an unambiguous quantity for describing the spin 
transport. We also address the addition dissipative Fermi 
surface contribution to spin polarization, which may be 
relevant to experimental measurement. 9 ' 10 

The remainder of the paper is organized as follows. 
In Sec. II, we numerically compute the SHC through 
the Kubo formula in a tight-binding model of the 2DHG 
system at different sample sizes and disorder strengths, 
and perform finite-size scaling analysis. In Sees. Ill A, 
IIIB and IIIC, we perform a Laughlin's gauge experi- 
ment numerically to determine the intrinsic spin trans- 



2 



fer/accumulation due to the adiabatic insertion of a mag- 
netic flux quantum. And in Sec. HID, we compare the 
results with a system in a perpendicular external mag- 
netic field and demonstrate the relationship between the 
intrinsic spin transfer (accumulation) and the existence 
of edge states. Finally, a summary is given in Sec. V. 



II. SPIN HALL CONDUCTANCE 

A. Hamiltonian 

We start with the 2DHG Hamiltonian proposed by 
Bernevig and Zhang 12 



For the convenience of the following numerical study, 
we further convert this continuum Hamiltonian into a 
tight-binding version on the square lattice. This may 
be realized by making the replacement k v — > sin k v and 
k 2 — > 2(1 — cos/c„). Two models are apparently equiva- 
lent near the band bottom where k v — > 0. The resulting 
Hamiltonian reads 

(ij) iv 

+ Vr j^y^^ c \Sy c i+% ~ i ^2 c iSx c *+y + H.c^j 
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FIG. 1: Band structure for a clean 2DHG system determined 
by Eq.( 2) (with t = 1). The HH and LH denote the heavy 
hole band and light hole band, respectively. 



The parameters r y 1 and -f 2 m GaAs are given by r y 1 — 6.92 
and 7 2 = 2.1, 19 so we have t ~ 1.45eV and V L ~ 0.345i. 
We shall choose Vr = 0.02i for the Rashba-type coupling. 
By assuming that the 2DHG is confined in a well of the 
thickness a ~ 8.7 nm, one approximately has (k z ) = 

and (k 2 z ) ~ {^j^j , where a$ ~ 0.565 nm for GaAs. The 

corresponding energy spectrum for a clean sample (e, = 
0) is shown in Fig. 1. Due to the confinement in the z 
direction, a gap with 0.029i <~ Q.i)A2eV opens up between 
the heavy hole (HH) band and the light hole (LH) band 
at r point (see Fig. 1). For simplicity, we shall set t = 1 
in the rest of paper. Finally, e» in Eq.(2) accounts for the 
on-site and spin-independent disorder strength, which is 
randomly distributed within [—W/2,W/2]. 



B. Spin Hall conductance 



(2) 



where the electron annihilation operator Cj has four 
components characterized by the "spin" index S z — 
|, i, — |, — |, respectively, and i + v (z> = x,y) denote the 
nearest- neighbors of site i. Here t is the nearest- neighbor 
hopping integral defined as 



2ma, 



2(7i + o7 2 ), 



(3) 



where m is the mass of the electron and ao is the lat- 
tice constant. Vl and Vr represents the strengths of the 
Luttinger-type and Rashba-type spin-orbital couplings, 
respectively: 



Vr = 



Vr = 



72 



7i + §72 



t 



7i 



§72 



The SHC is defined by the Kubo formula 18 



\ E n <E f <E„ 



{n\j^\rn){m\j y \n) 



(E m — E n ) 



(4) 



where N is the total number of lattice sites, Ef denotes 
the Fermi energy, E m (E n ) is the eigen-energy, and (...) 
is averaged over all disorder configurations. The charge 
current operator and spin current operator j^ spin are 
defined as 



1 

r 



where the velocity operator v is the conjugate operator of 
the position operator R = ^2 ir7 ^72^ (n icr is the electron 
number operator at site i with spin index a), is defined 
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FIG. 2: as£r vs Ef at VF = 0. The inset shows the whole 
Fermi energy region and the main panel shows the details 
around the band bottom. The solid line is at 1000 x 1000 
lattice, the dashed line is at 24 x 24 lattice, and the open 
circles are at 24 x 24 lattice with averaged over 1000 BCs. 
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by the standard relation 



= jr[H,R»] 



(5) 



For a pure sample (ej = 0), the Hamiltonian can be 
easily diagonalized in momentum space, and we can de- 
termine ctsh vs Ef for a very large lattice (N — L 2 with 
L = 1000) as shown in the inset of Fig. 2. Since the 
tight-binding model is a good approximation of the orig- 
inal model only when Ef is near the band bottom, we 
shall focus on the regime of < Ef < 0.6 as represented 
in the main panel of Fig. 2 by the solid curve. One can 
see that <tsh in this region is a little bit larger than 10 
e/87r and very flat as a function of Ef except for that 
within the gap (~ 0.029) between the HH and LH bands. 

Because the momentum is no longer a good quantum 
number in the presence of disorder, Eq.(2) can not be 
diagonalized in the momentum space when =/= 0. In 
the latter case, the maximal lattice size reachable by the 
numerical calculation is usually much smaller than L = 
1000 (around L = 20 — 30). To check what happens at 
smaller lattices, we compute asn on an L = 24 lattice in 
the pure system with a periodic boundary condition (BC) 
and the result is illustrated by the dotted curve in Fig. 2. 
The corresponding asH exhibits very large fluctuations 
whose magnitude can be several times larger than the 
converged value at L = 1000. To remove such finite-size 
fluctuations, we introduce an average over different BCs. 
A general (twisted) BC is given by 



ijj(x + L,y) 
ip(x,y + L) 



^(x,y) 
"ip(x,y) 



(6) 



where # M 's are defined within [0, 2ir] with 6^ = or 2-7T 
(p = x, y) corresponding to the periodic BC. Since the 



FIG. 3: a S H vs. Ej at W = 0.2, 1, and 3. The solid curves 
are for L = 16, and the dashed ones for L — 24. (a) is for 
the 2DHG system and (b) is for the spin-3/2 Rashba model 
by setting Vl = 0. 



SHC should not be sensitive to the BCs in the thermo- 
dynamic limit, in principle, ush in Eq. (4) can be al- 
ways defined as averaged over both disorder and twisted 
BCs in such a limit. For a clean system, the averaged 
<7sh at L = 24 over 1000 randomly generated BCs is 
shown by the open circles in Fig. 2. One can clearly 
see that the large fluctuations in asH at L — 24 (dotted 
curve) have been significantly reduced after the BC av- 
eraging, and the result (open circles) coincides with asn 
obtained at L = 1000 (solid line) very well. It indicates 
that the finite-size effect at smaller size (L = 24) can 
be properly removed by the BC average. As a conse- 
quence, the BC averaged SHC behaves much smoother 
than at a fixed BC, and is thus more suitable for a finite- 
size scaling analysis. We note that the BC average is also 
required by the gauge- invariant condition for SHC in a fi- 
nite size system. When we apply an electric field through 
a time-dependent vector potential, it acts as a gener- 
alized boundary phase evolving with time (see below). 
Thus the SHC averaged over time is gauge-invariant and 
equivalent to the BC averaged SHC. 6 Such a numerical 
method has been previously used in the study of the 2D 
Rashba model and 3D Luttinger model. 6,7 Similar bound- 
ary phase averaged charge Hall conductance can be also 
related to a topological invariant Chern number. 18 

Now we use the above method to study the sample- 
size dependence of asn at W — 0.2 for weak disorder 
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FIG. 4: W and sample-size dependence of <jsh with aver- 
aged within the interval [0.3,0.4]. (a) gsh vs. W. The open 
squares are for L = 12, close squares for L = 16, open circles 
for L = 20, and close circles for L — 24. (b) ctsh vs L. The 
lines from top to bottom are for W = 0.05, 0.2, 0.6, 1, 1.5, 2, 
2.5, 3, and 5, respectively. 



strength; W = 1 for intermediate disorder; and W = 3 
for strong disorder. The calculated ctsh is presented 
in Fig. 3 (a) for two sample sizes: L = 16 [averaged 
over 20000 configurations (solid line)] and L = 24 [aver- 
aged over 5000 configurations (dashed line)]. The size- 
dependence of asH is very weak for all three cases at 
different -E/'s. For comparison, we have also done a sim- 
ilar calculation for a "spin-3/2 (linear) Rashba model", 
which corresponds to Eq. (2) at Vl = and Vr = 0.5. 
SHC at W = 0.2, 1, and 3, with L = 16 and 24, respec- 
tively, arc shown in Fig. 3(b), which clearly decrease with 
the increase of the sample length, similar to the result in 
the spin- 1/2 Rashba model. 6 

To further investigate the effects of the lattice size and 
disorder, we calculate the averaged asH within a chosen 
energy interval [0.3, 0.4] at various disorder strengths and 
sample sizes as shown in Fig. 4. <jsh is almost indepen- 
dent of sample sizes and reduces monotonically with in- 
creasing W, which is different from the spin-1/2 Rashba 
model but similar to the behavior in the 3D Luttinger 
model. 7 In the weak disorder limit, ctsh approaches to 
the value slightly larger than 10 e/87r, which is quite close 
to the value (~ 9e/87r) obtained 20 in an analytic calcu- 
lation in a different model for 2DHG system. Thus we 
establish that SHC is finite extrapolating to thermody- 
namic limit until the disorder strength reaches a critical 
value W c = 5.0. 



III. INTRINSIC SPIN TRANSFER AND 
ACCUMULATION 

In the last section, we have shown the numerical ev- 
idence that a finite <tsh can survive in the presence of 
disorder with the sample size being extrapolated to the 
thermodynamic limit. This is in sharp contrast to the 
2D Rashba model in which asH vanishes in the thermo- 
dynamic limit with the turn on of very weak disorder 
strength. 6 Such a distinct behavior in osh is consistent 



FIG. 5: (a) The geometry of the sample in the gauge exper- 
iment. It is like a belt enclosing a flux in the center, (b) 
The implementation of (a) in the numerical calculation. The 
dashed line with a phase factor e l * along the x axis means a 
twist BC, which represents the flux in (a). 

with the perturbation theories in which the vertex cor- 
rection also results in opposite conclusions on the fate of 
gsh in two models, as mentioned in the Introduction. 

In general, a finite conductance should lead to an ac- 
cumulation of current carriers at the edges (open bound- 
aries) of the sample according to the current conservation 
law. However, in the SOC system the spin current con- 
servation law does not exist as the spin is not conserved. 
Thus the relation between the SHC and the spin accu- 
mulation is not straightforwardly present. 

In this section, we will try to probe the spin transfer 
or accumulation effect directly by performing numerically 
a Laughlin's gauge experiment, which was first used in 
the IQHE system 21 ' 22 to explain why an integer num- 
ber of charges can be transported across the sample in 
the transverse direction when a longitudinal electric field 
is applied via an adiabatic flux insertion process. This 
method has been recently generalized to the SOC sys- 
tem by Sheng et al. 6 in the study of the SHE in the 2D 
Rashba model. 



A. Gauge experiment and level crossing 

The numerical gauge experiment will be performed on 
a belt-shaped sample with a magnetic flux $ adiabati- 
cally inserted at the center of the ribbon [see Fig. 5(a)]. 
Such a belt-shaped geometry with the magnetic flux can 
be implemented in a square lattice sample shown in Fig. 
5(b) by imposing a twisted BC along the x direction and 
an open BC along the y direction. Namely, with a gauge 
transformation, one may impose the flux <j> by changing 
the BC along x direction from the periodic BC to a twist 
BC with the phase twist 6 X = $ indicated by the dashed 
line and the phase factor e i$ in Fig. 5(b). 

Numerically the energy spectrum of single electron 
states at a fix $ can be calculated in the geometry of 
Fig. 5(b). Fig. 6 shows such calculated energy spec- 
tra as a function of $ in the pure system. Note that 
the system is equivalent at $ = and <j> = 2tt as the 
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FIG. 6: The evolution of the single electron state energy 
with $ for the 2DHG on a 32 x 32 clean sample (W = 0). 



Hamiltonian is periodic = 0) = H(& — 2tt), and 

the energy spectrum is symmetric between < <f> < tt 
and 7r < $ < 2n due to a time-reversal symmetry, such 
that only the first half is plotted in Fig. 6. 

One may imagine changing the flux <I> adiabatically 
from to 2n to follow the corresponding change of the 
system. Generally there arc two possibilities for the evo- 
lution of the ground state with the electrons filled be- 
low the Fermi energy Ef. If each single electron state 
always anti-crosses with other states (level repulsion), it 
will eventually return to the original state after the inser- 
tion of one flux quantum of <I> = 2tt sufficiently slow (i.e., 
adiabatically such that there is no quantum transition to 
other levels during this operation); On the other hand, 
if a level crossing occurs for single electron states during 
the adiabatical increase of $, then the original ground 
state generally will evolve into a different state through 
the change of the single state occupation number, even 
though the single electron energy spectrum must remain 
the same after the insertion of one flux quantum. 

A slow increase of the flux $ is equivalent to applying 
a weak electric field along the x direction: 



frr = - 



1 d<S> 
L x dt ' 



(7) 



Then the intrinsic charge/spin Hall effect due to an elec- 
tric field in belt-shaped system will lead to the transfer 
of charge/spin along the y-axis during the flux insertion. 




FIG. 7: The miniband structure of a 32 x 32 clean sample. 
Two arrows represent the evolution of two states by tuning 
the parameter which start from the beginning at $ = 
and arrive at the end of the arrows at $ = 2ir (see the text 
for details). The dashed line indicates the Fermi energy. 



If there is to be charge/spin accumulation at the open 
edges [Fig. 5(a)], it means that the ground state should 
not return to the original one after the adiabatic flux 
insertion, which is the essence of Laughlin's gauge exper- 
iment for the IQHE system. Therefore, the level crossing 
of single electron states serves as a necessary condition 
for the existence of intrinsic spin transfer/accumulation 
for a SHE system. As clearly shown in Fig. 6, which 
is determined on an L = 32 lattice, the single electron 
energy levels do simply cross in the pure 2DHG system. 
The resulting net spin transfer will be analyzed in the 
following section. 



B. Spin Transfer at W = 

In this section, we calculate the spin transfer induced 
by tuning the flux $ from to 2ir in the pure case 
(ei = 0). Starting from the periodic BC ($ = 0), the 
momentum k x along x direction is a good quantum num- 
ber. And due to the open boundary along the y direction, 
the electron states form a mini-band structure which is 
shown in Fig. 7 and can be denoted as \k x ,n), where 
n is the index of the mini-bands. If we choose an arbi- 
trary Fermi energy Ef, the ground state corresponds to 
that all the single particle states below Ef are occupied, 
whereas unoccupied above Ef. So the spin density along 
y direction (labelled by y = j y ) at the initial state can 
be defined by 



Sl(y) 



(k x ,n\^2c+S z Cj\k x ,n). (8) 



E 

E kx , n <E f 



When one inserts the flux adiabatically from zero to $, 
the state \k x ,n) will evolve to \P X = k x + Q/L x ,n), For 
k x > and < 0, the corresponding energies of the states 
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will increase and decrease, respectively, which is illus- 
trated by arrows in Fig. 7. Thus one finds a simple level 
crossing for two levels close to each other with opposite 
signs of k x , which represents the level crossing shown in 
Fig. 6 for pure system due to no mixing term between 
different states in the Hamiltonian. After the insertion 
of one flux quantum, the state \k x ,n) evolves into a new 
state \k x + 2n/L x ,n) at $ = 2tt. Then the spin density 
at the final state is given by 

s f Av)= E <fc* + f^E4^#. + §U}. 

Ek x , n <E f x j x 

(9) 

The resulting change in the spin-density due to threading 
a flux quantum is AS z (y) = Sl(y) — S l z (y). 

AS z (y) in a pure 2DHG with 500 x 1 00 lattice is shown 
in the inset of Fig. 8(a), where the Fermi energy is 
Ef = 0.4. We recall that due to the time-reversal invari- 
ance, S\{y) — everywhere. After adding the 2tt flux 
quantum, the spin density shows some strong peaks at 
the boundaries whereas with a much weaker magnitude 
in the bulk. We show the spin-wave like modulation in 
AS z (y) near one edge of the belt-shaped sample at var- 
ious lattice sizes in the main panel of Fig. 8(a). Here 
AS z (y) is almost sample-size independent at the edge, 
whereas it reduces monotonically with the sample size in 
the bulk. So in the thermodynamic limit, there should be 
only the spin density modulation present near the edge. 
Similar results with smaller magnitude are also obtained 
for an electron spin-f/2 Rashba model as shown in Fig. 
8(b) for comparison. Such a spin density modulation 
does support a net spin transfer along the y direction, 
accompanying an electric field generated along the x di- 
rection due to the inserting flux according to Eq. (7), 
but the magnitude of the total (integrated) AS Z ~ 0.2h 
at one edge is much less than one spin quantum, which 
is in sharp contrast to the large bulk SHC (the calcu- 
lated <jsh ~ lle/87r, that suggests a total spin transfer 
of 5.5(7i/2) spin quanta). This large discrepancy between 
two approaches indicates that the SHC is no longer an 
unambiguous quantity for measuring the spin transport 
even in the pure case. In the following we further explore 
the disorder case. 
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y 

FIG. 8: The spin density changes between the state at the 
beginning and the end of the procedure of increasing the flux 
from to 2ir in the pure samples. The main panel shows the 
data around one edge for various lattice sizes, while the inset 
shows the whole result on a 500 x 100 lattice. The Fermi 
energy is fixed at Ef = 0.4. (a) is for the 2DHG; The solid 
lines with open squares are for 1000 x 100 lattice, the dashed 
lines with close squares are for 500 x 500 lattice, and the dotted 
lines with open circles are for 1000 x 1000 lattice, (b) is for 
the 2D Rashba model; The solid lines with open squares are 
for 100 x 100 lattice, the dashed lines with close squares are 
for 500 x 100 lattice, and the dotted lines with open circles 
are for 1000 x 1000 lattice. 



C. Disorder effect 

Now we consider the disorders effect. Fig. 9 shows the 
results for a chosen disorder configuration with W = 0.05 
in a 32 x 32 lattice. A careful examination of the en- 
ergy levels reveals that each eigenstate goes up and down, 
making several large angle turns due to backward scat- 
tering. These energy levels never cross with each other, 
except for at <!> = and it, where two levels become 
exactly degenerate (Kramers degeneracy). Then if one 
follows any Kramers degenerated pair of states starting 
from $ = to $ = 2ir, one will always go back to exactly 
the same pair of states. Namely, after the adiabatic inser- 



tion of one flux quantum, all states evolve exactly back 
to the starting states. 

For comparison with the W = case, in Fig. 10(a) we 
combine an enlarged plot of two adjacent energy levels in 
Fig. 10 together with the pure case, where the solid lines 
are for the disordered case and the dashed lines are for 
the clean case. In the pure case, an electron at the state 
\h) (K2)) will finally evolve to (I/2)) with increasing 
flux due to the level crossing. However, in contrast to 
the level crossing in the pure system, two energy levels, 
|«2), in the disordered system generally show a level 
repulsion as one increases the flux <!>, to evolve into the 
final states of I/2} and respectively [see Fig. 10(a)]. 
The level repulsion in the disordered case represents a 
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FIG. 9: The evolution of the single electron states with <E> 
for the 2DHG on 32 x 32 lattice for a disordered sample with 
W = 0.05. 




FIG. 10: (a) The comparison of the evolutions of two levels 
for clean system (dashed line) and disordered system with 
W = 0.05 (dashed line). The levels are crossing at clean 
system and anti-crossing at disorder system. The and 
\i2) are the initial states while the and I/2) are the final 
states, see the text for details, (b) The process of a state 
scattered to I/2} by the disorder. Such a process is leading to 
the anti-crossing as shown in (a). 



scattering process shown in Fig. 10(b). 

We have done similar calculations at different lattice 
sizes, such as 8x8, 16x16, and 24 x 24. The results 
are all similar except that the "gaps" characterizing the 
level repulsion decrease with the lattice size, which are 
roughly proportional to 1 /N. According to the discussion 
in Sec. Ill A, the adiabatic insertion of the flux <E> will 
not generate true spin transfer or accumulation as the 



FIG. 11: The evolution of the single electron states with $ in 
the system with an external magnetic field on a 32 x 32 lattice. 
The strength of the magnetic field is 7r/8 per plaquette. 



ground state simply evolves back after a flux quantum 
insertion. In fact, we have checked various weak disorder 
strengths and always found the same generic result. 

Therefore, we conclude that there is no intrinsic spin 
Hall accumulation (spin transfer) for the 2DHG system 
in the disorder case. This conclusion is in sharp con- 
trast to the calculated SHC at small W in Sec. II, which 
remains finite and close to the value of the pure case 
when the disorder is weak. According to the discussion 
at the beginning of Sec. Ill, these two approaches, for a 
bulk system on a torus or as an open system with two 
edges, are not necessarily in accordance with each other, 
because the spin current continuity equation is broken. 
For example, a backward scattering of an electron by an 
impurity shown by Fig. 10(b) can simultaneously change 
the sign of its momentum as well as its spin direction due 
to the SOC. Then one can find that the spin current is 
still additive (same sign) in this scattering process, which 
contributes to the SHC accordingly, but there should be 
no true spin transfer taking place because of the back- 
ward scattering of the spin carrier - electron which moves 
back in direction after the scattering. 



D. Edge state 

To further understand the condition for the presence of 
level crossing in the disordered case, we apply an exter- 
nal magnetic field perpendicular to the 2D sample (with 
a flux strength 7r/8 per plaquette) and repeat the above 
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calculation. We find that most of the "gaps" still behave 
as 1/N , similar to those in the non-magnctic-ficld case. 
These are bulk states. But we can also identify a few 
levels whose level-repulsion "gaps" reduce exponentially 
with TV with the increase of the sample size. For exam- 
ple, such kind of "gap" at 32 x 32 lattice size is about 
1/30 of the one at 16 x 16 lattice size, an exponential 
dependence on system length. We show some of these 
"level crossings" in Fig. 11 with the same parameters 
as in Fig. 9. The "gaps" in Fig. 11 are much smaller 
than the ones in Fig. 9 and hard to distinguish with the 
naked eye. The exponential decaying form of the gaps 
indicate the existence of level crossings in the thermo- 
dynamic limit in the presence of a perpendicular mag- 
netic field, which corresponds to an IQHE. Then a net 
charge (spin) is transported from one edge to the other 
correspondingly. 6 

We have also systematically compared the spatial dis- 
tribution of electron density p i = (cjc,) of the states in 
the 2DHG system with and without the presence of a 
perpendicular magnetic field. It is found that the level 
crossing only occurs for the edge states. Physically this 
is because in the thermodynamic limit the distance be- 
tween two edge states in the opposite sides of the sample 
goes to infinity which makes the scattering between them 
vanishes exponentially Thus the absence of level cross- 
ing (and the intrinsic spin Hall transfer) in the 2DHG 
system can be attributed to the fact that there do not 
exist any edge states in such a system in the absence of 
perpendicular external magnetic fields. 

IV. CONCLUSIONS AND DISCUSSIONS 

In the present work, we have numerically calculated the 
SHC through the Kubo formula at different lattice sizes 
and impurity strengths. A reasonable finite-size scaling 
analysis shows that the bulk SHC remains finite in the 
thermodynamic limit in the disordered 2DHG up to a 
critical disorder strength W c — 5. This result is con- 
sistent with the perturbative calculation based on the 
vertex correction. 12 

However, we have also found the strong evidence that 
the net spin Hall accumulation or spin transfer disappears 
in the disordered case based on a numerically performed 
Laughlin's gauge experiment. We have interpreted the 
conflicting results of the SHC and direct spin accumula- 
tion calculations as due to the fact that the spin is not 
conserved in an SOC system and there is no longer a 
conservation law to govern the spin current. Therefore 
the SHC is no more a unique quantity to characterize the 
spin transport, where electron spin can relax and vanish 
in the bulk without being transported across the sample 
by a transverse electric field. 

By comparing with the case in the presence of an ap- 
plied perpendicular magnetic field, which is in the IQHE 
regime, we have shown that the absence of an intrinsic 
spin transfer (accumulation) can be related to that there 



are no edge states in the 2DHG system, in contrast to 
the IQHE case. 

We further address the validity of the above adiabatic 
argument in the thermodynamic limit. It is noted that 
since an intrinsic spin transfer (accumulation) is always 
absent in a finite size system according to the above ar- 
gument, it is very hard to imagine that such an effect 
could be restored as an intrinsic one by the Landau- 
Zener tunnelling 8 between two anticrossing levels in the 
thermodynamic limit. Of course, the Landau-Zener tun- 
nelling between the level repulsion gaps in the large sam- 
ple size limit (the level repulsion gap vanishes in a 1/N 
fashion, with the same rate as the average level spacing 
vanishing in this limit) can still contribute to dissipa- 
tive transport currents, and similar to the conventional 
charge transport, one may expect a dissipative spin trans- 
port term appearing, besides the null intrinsic spin trans- 
fer. Here the dissipative spin transport is qualitatively 
different from the true level crossing with edge states. 
At a finite W, in the longitudinal channel, the so-called 
"Thouless conductance" term can be obtained based on 
the sensitivity of the energy level to the threading flux, 
giving rise to a finite a xx , 23 which is subject to the de- 
tailed scattering mechanism as contributed by the levels 
near the Fermi energy. In the present SOC system, a 
dissipative spin Hall transfer related to this longitudinal 
conductance is indeed possible, which similarly involves 
only the states near the Fermi energy. We note that such 
a dissipative term may also be important for the spin- 
polarization measured experimentally 9 ' 10 . 

The vanishing spin transfer revealed through Laugh- 
lin's Gcdanken gauge experiment for the present 2DHG 
and 2D electron Rashba model 6 under an adiabatic pro- 
cedure may be a generic behavior for a metallic system. 24 
It would be interesting to directly probe this property for 
more SOC systems, including the 3D Luttinger model, al- 
though numerically it is more challenging for the latter as 
it is much harder to do the finite-size scaling for the level 
repulsion gaps in a 3D system. Finally, it is interesting to 
note that several recent papers have studied topological 
SHE in insulating electron systems with SOC. 25 In these 
models, it becomes clear that bulk intrinsic SHC can re- 
sult in spin transfer and accumulation, with the presence 
of spin-polarized edge states, which is associated with a 
topological invariant. This is in agreement with our gen- 
eral conclusion that the intrinsic spin transfer relies on 
the existence of current-carrying edge state. 
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